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Chapter  I  Exterior  Algebra 


Introduction 

This  report  will  bring  together  several  mathematical  topics  which 
seem  to  be  relevant  to  the  treatment  of  nonlinear  networks.  The  main 
goal  will  be  to  discuss  some  stability  theorems  recently  obtained  by 
Brayton  and  Moser  [l].  Before  reaching  this  goal,  we  shall  present 
some  exterior,  or  Grassman  algebra,  some  algebraic  topology,  some 
theorems  about  electrical  networks,  and  a  few  remarks  about  Liapunov 
stability.  Of  course,  the  final  results  could  be  obtained  more  directly 
without  the  excursions  into  exterior  algebra  or  algebraic  topology. 
However,  since  these  topics  do  give  some  additional  insight,  it  was 
thought  worthwhile  to  include  them.  Because  there  are  many  places 
where  a  rigorous  and  more  complete  treatment  may  be  found,  the  approach 
throughout  has  been  intentionally  kept  to  a  verv  elementary  and  a  very 
intuitive  level. 

Linear  Spaces 

We  shall  begin  with  a  discussion  of  the  algebra  of  three-dimensional 
vectors.  Mathematicians  have  learned  that  it  is  more  important  to  know 
how  things  behave,  than  to  know  what  they  are.  One  well-known  fact 
about  the  behavior  of  vectors  is  that  they  can  be  added,  that  is,  two 
vectors  a  and  b  can  be  combined  to  form  a  third  vector  which  is 
called  their  sum  and  which  is  written  a  +  b.  The  operation  by  which  the 
sum  is  formed  has  the  following  properties: 

The  associative  property,  that  is,  a  +  (b  +  c)  =  (a  +  b)  +  c ; 


1. 


2 


2.  a  neutral  element,  written  0,  exists  such  that  for  all 
vectors  a,  we  have  a+0  =  0  +  a  =  a; 

3-  every  vector  a  has  an  inverse  element  -  a  such  that 
a  +  (-  a)  =  (-  a)  +  a  =  0; 

4.  the  commutative  property,  that  is,  a  +  b  —  b  +  a. 

These  properties  are  evident  from  the  geometrical  definition  of 
addition  by  the  Parallelogram  Law  which  is  illustrated  in  Figure  1. 


Any  set  of  elements  for  which  an  operation  such  as  +  can  be  carried 
out,  and  for  which  the  properties  1,2,3,  are  satisfied  is  called  a 
g roup .  If  Property  4  is  also  satisfied,  the  group  is  said  to  be  com¬ 
mutative. 


Besides  addition,  there  is  another  operation  possible  with  vectors, 
namely  multiplication  with  a  scalar.  Let  us  denote  scalars  by  Greek 
letters ;  then  the  operation  of  multiplying  by  a  scalar  has  the  following 
properties : 


5. 

a(a  +  b)  =  aa  +  ab; 

6.  (a  +  p)a  =  aa  +  pa; 

7. 

a(pa)  =  (ap)aj 

8  .  la  =  a . 
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Any  set  of  elements  for  which  an  operation  having  the  properties 
l,2,3j  and  4  can  be  defined,  and  for  which  an  operation  having  the 
properties  5,6,7,  and  8  can  be  defined  is  called  a  linear  space ,  or  a 
linear  vector  space.  The  elements  a^,a2,...,an  in  a  linear  space  are 
said  to  be  linearly  independ ent  if  a  relation  such  as 


a. a1  +  a„a„  +•••+  a  a  =  0 
11  2  2  n  n 


implies  that  =  • • •  =  an  =  0 .  If  the  elements 

are  not  linearly  independent,  they  are  said  to  be  linearly  dependent. 


The  geometrical  significance  of  linear  independence  is  clear.  One 
vector  is  linearly  independent  if,  and  only  if,  it  is  not  the  zero  vector. 
Two  vectors  are  linearly  independent  if,  and  only  if,  they  do  not  lie  in 
the  same  line.  Three  vectors  are  linearly  independent  if,  and  only  if, 
they  do  not  lie  in  the  same  plane. 


The  vectors  a^,a^, - • • >an  are  said  to  form  a  basis  of  the  linear 
space  if  every  element  of  the  space  can  be  expressed  uniquely  as  a  linear 
combination  of  the  basis,  that  is,  given  any  element  x  in  the  linear 
space,  there  exists  a  unique  set  of  scalars  ci^q^,  •  •  •  >an  such  that 


x  =  a1&1  +  a2a2  +• 


•+  a  a  . 
n  n 


The  number  of  vectors  in  the  basis  is  called  the  dimension  of  the  space. 

The  space  of  geometry  and  physics  obviously  has  dimension  three  because 

every  vector  can  be  expressed  uniquely  as  a  linear  combination  of  the 
— ►  — ►  — ► 

fundamental  i , j , k  vectors . 
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Bivectors 

Besides  the  operations  of  addition  and  of  multiplying  by  a  scalar, 
there  are  other  vector  operations  possible,  namely  the  operations  of 
multiplying  two  vectors  together  to  obtain  the  dot  product  and  the  cross 
product.  In  this  section  we  shall  introduce  a  third  type  of  product, 
the  exterior  or  wedge  product.  The  wedge  product  is  very  closely  related 
to  the  vector  product  in  a  three-dimensional  space.  However,  in  higher 
dimensional  spaces,  the  vector  product  does  not  exist  but  the  wedge  product 
always  exists.  Also,  the  wedge  product  has  the  associative  property 
which  the  vector  product  does  not. 

We  shall  introduce  the  wedge  product  geometrically.  A  vector  has 
the  geometrical  significance  of  a  directed  line  segment  which  has  a  magnitude, 
its  length,  and  a  direction.  Consider  a  parallelogram  (see  Figure  2)  having 
the  vectors  a  and  b  as  sides.  This  parallelogram  has  a  magnitude, 


its  area,  but  we  shall  see  that  it  is  useful  to  assign  to  the  parallelo¬ 
gram  also  an  orientation,  that  is,  a  sense  or  direction  in  which  the  sides 
are  traversed.  To  see  the  need  for  an  orientation,  consider  a  triangle 
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(see  Figure  3)  whose  area  we  shall  denote  by  ABC.  Let  0  be  a  point 


C 


Figure  3 


inside  the  triangle,  then 

ABC  =  0A3  +  OBC"  +  OCA . 


(1) 


However,  if  the  point  0  is  outside  the  triangle  as  in  Figure  4, 


then  the  equation  (1)  is  no  longer  true.  Instead,  we  have 

ABC  =  OAeF  -  03C  +  OCA.  (2) 


It  would  certainly  be  convenient  to  be  able  to  use  equation  (l)  no 
matter  where  0  is  placed.  This  can  be  done  if  we  agree  that  an  area 
is  to  be  considered  positive  or  negative  according  as  its  vertices  are 
described  in  a  counter-clockwise  or  a  clockwise  order.  For  example, 

ABC  will  denote  the  positive  number  equal  to  the  magnitude  of  the  area 
of  the  triangle  in  Figure  3  or  Figure  45  whereas  ACB  will  denote  the 
negative  number  whose  absolute  value  is  equal  to  the  magnitude  of  the 
area  of  the  triangle.  With  this  interpretation,  equation  (l)  is  valid 
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for  both  Figures  3  and  4,  because  in  Figure  3  OBC  describes  the  triangle 
in  a  counter-clockwise  sense  and  thus  OBC  represents  a  positive  number, 
whereas  in  Figure  4  OBC  describes  the  triangle  in  a  clockwise  sense 
and  thus  OBC  represents  a  negative  number.  Henceforth,  we  shall  con¬ 
sider  all  areas  to  be  oriented. 

Just  as  two  vectors  are  considered  equal  if  they  have  the  same 
length  and  are  parallel  in  the  same  direction,  so  we  shall  consider  two 
parallelograms  equal  if  they  have  equal  areas,  the  same  orientation,  and 
lie  in  the  same  or  parallel  planes.  We  shall  say  these  parallelograms 
represent  a  bivector.  We  shall  denote  by  aAb  a  bivector  such  as  the 
one  represented  by  the  parallelogram  in  Figure  2,  if  the  orientation  of 
the  parallelogram  is  such  that  the  vector  a  must  be  rotated  in  the 
counter-clockwise  sense  to  coincide  with  the  vector  b.  However,  if  this 
parallelogram  is  oriented  in  the  clockwise  sense,  we  denote  the  bivector 
by  bAa.  Because  the  only  difference  between  the  bi vectors  is  in  their 
orientation  or,  what  is  equivalent,  in  the  sign  attached  to  the  area,  we 
shall  write 

aAb  =  -  bAa.  (3) 

Thus,  the  multiplication  indicated  in  the  wedge  product  is  anti -commutative 
just  as  is  the  multiplication  in  the  cross  product.  Notice  that  since 
aAa  =  -  aAa,  we  must  have  2 aAa  =  0,  or  aAe  =0.  We  conclude  that 
the  wedge  product  of  a  vector  with  itself  is  always  zero. 

It  should  be  emphasized  that  a  wedge  product  of  two  vectors  is  not 
a  vector  but  it  is  a  bivector  which  is  represented  by  a  parallelogram 
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with  a  definite  orientation.  We  can  introduce  a  procedure  for  adding 
two  bi vectors  together.  This  addition  procedure  is  indicated  in  Figure  5. 


Figure  5 

We  find  that  this  addition  has  the  properties  1, 2,3,4  of  the  preceding 
section. 

We  can  also  define  the  operation  of  multiplying  a  bivector  by  a 
scalar.  The  product  of  the  bivector  aAb  by  the  scalar  a,  written 
caAb,  is  the  bivector  represented  by  a  parallelogram  whose  area  is 
| a |  times  the  area  of  the  parallelogram  in  Figure  2  representing  aAb, 
whose  orientation  is  the  same  or  opposite  to  that  in  Figure  2  according 
as  a  is  greater  or  less  than  zero,  and  whose  plane  is  the  same  or  parallel 
to  that  of  Figure  2.  It  is  easy  to  see  that  this  operation  of  multiplying 
a  bivector  by  a  scalar  has  properties  5, 6, 7, 8  of  the  preceding  section. 
Consequently,  the  set  of  bivectors  forms  a  linear  space.  We  shall  soon 
show  that  this  space  is  three-dimensional. 

Let  us  illustrate  the  procedure  for  finding  a  wedge  product  by 
considering  a  simple  example.  Suppose  a  =  a^i  +  a^j  +  a^k  and 
b  =  P-jl  +  P2j  +  P^k,  then 


aAb  =  a1P2  iAj  +  a ^  iAk  +  jAi  +  jAk  +  a ^  kAi  +  a^2  kAj .  (4) 
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Notice  we  have  used  the  fact  that  the  wedge  product  of  a  vector  with 
itself  is  zero.  The  result  can  be  still  further  simplified  by  using  the 
anti -commutativity  (3)  of  the  product.  We  have  jAi  =  -  iAj,  kAj  =  -  jAk, 
iAk  =  -  kAi,  and  therefore  the  product  in  (4)  simplifies  to 

^■Ab  ( ~  a2^l^  ^ ^p^3  —  “3^2^  J ^3^1  "  ^1^3)  kAi.  (5) 

Since  every  vector  can  be  written  as  a  linear  combination  of  i,j,  and 
k,  it  is  clear  that  every  bivector  will  be  a  linear  combination  of  the 
three  bivectors  iAj,  jAk,  and  kAi.  This  shows  that  the  space  of  bi vectors 
is  three-dimensional. 

The  equation  (5)  suggests  the  connection  between  the  wedge  product 
and  the  vector  product.  Notice  that  each  of  the  fundamental  bi  vectors 
iAj,  jAk,  and  kAi,  has  a  unique  normal,  namely,  the  vectors  k,i,  and  j, 
respectively.  If  we  replace  each  bivector  in  (5)  by  its  normal,  we  get 
the  vector 

~  a3|32^1  +  (a3^i  ~  “1^3) J  +  (“^^2  "  a2Pl^k’ 
which  is  the  vector  product  of  a  and  b. 

Determinants 

The  procedure  for  forming  wedge  products  can  be  extended  to  form  the 
wedge  product  of  three  vectors.  This  triple  wedge  product,  which  we  call 
a  tri vector,  is  represented  geometrically  by  the  oriented  volume  of  the 


parallelepiped  generated  by  the  three  vectors  (see  Figure  6) .  The  volume 
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is  oriented  by  assigning  an  order  in  which  the  vectors  a,b,  and  c  are 


Figure  6 


described.  If  the  order  is  such  that  the  vectors  form  a  right-handed 
set,  the  volume  is  considered  positive.  If  the  vectors  form  a  left- 
handed  set,  the  volume  is  considered  negative. 

Let  us  illustrate  the  procedure  for  forming  the  tri vector  by  con¬ 
sidering  an  example  in  which  the  vectors  a  and  b  are  the  vectors 
used  in  (4)  and  the  vector  c  =  Y^i  +  Ygj  +  Y^k*  Using  (5),  we  find 

that 

aAbAc  =  Y3(aip2  -  a^)  iAjAk  +  Y-^a,,^  -  a^z)  jAkAi  +  Y-^a^  -  a ^  )  kAiAj. 

(6) 

Because  of  (3) ,  we  see  that  kAiAj  =  -  iAkAj  =  iAjAk,  and  jAkAi  =  iAjAk, 
therefore  the  result  in  (6)  can  be  simplified  to  the  following: 

aAbAc  =  [Y1(a2P-3  -  Y2^°3^1  “  r^al^2  "  °2^1^  lAJAk. 

The  expression  in  the  bracket  is  obviously  the  determinant  of  the 
numbers  which  are  the  components  of  the  vectors  a,b,  and  c.  This  result 
justifies  the  geometrical  interpretation  of  the  trivector  as  the  volume 
of  the  parallelepiped  in  Figure  6. 
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A  wedge  product  for  any  two  vectors  with  n  components  can  be 
defined  by  a  slight  extension  of  the  method  used  to  define  the  wedge 
product  of  two  vectors  in  three-dimensional  space.  For  example,  if 

a  =  (2,3,-  1,5),  b  =  (0,1,-  2,4) 
are  two  4-vectors,  we  first  write  them  as 

a  =  2e1  +  3e2  -  e^  +  5e^,  b  =  e2  -  2e^  +  4e^ 

where  e^,e2,e^,e^  are  the  basic  4-vectors  defined  by  the  equations 

ej  —  (1,0, 0,0),  e2  —  (0, 1,0,0), 

e^  —  (0,0, 1,0),  -  (0,0, 0,1). 

We  then  form  the  wedge  product  using  the  anti -commutativity  rule  (3)  and 
the  fact  that  the  wedge  product  of  two  identical  vectors  is  zero.  We  get 

aAb  =  2  e^Ae2  -  4 e^Ae^  +8  e^Ae^  -  5  e2Ae^  +  7  e2Ae/  +  ^  e^Ae^. 

The  wedge  product  of  any  number  of  n-vectors  can  be  defined  by 
continued  multiplication.  For  example,  if 

c  =  el  -  e2  +  3e4, 

d  =  2e.  -  2e„  +  e,, 

X  3  4 

then  the  wedge  product  of  the  vectors  a,b,  and  c  is 

aAbAC  =  -  9  e]_Ae2Ae3  +  ^1  e^Ae^e^  -  6  e^Ae^Ae^  -  21  e2Ae3Ae4 
and  the  wedge  product  of  the  vectors  a,b,c,  and  d  is 
aAbACAd  =  75  e^AegAe^Ae^. 
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The  number  75  is  the  value  of  the  determinant  of  the  four  vectors  a,b,c, 
and  d,  that  is, 


This  connection  between  the  determinant  of  four  4-vectors  and  their 
wedge  product  suggests  the  following  theorem:  An  n  by_  n  determinant 
may  be  evaluated  by  forming  the  wedge  product  of  the  n  column  vectors . 
Of  course,  the  determinant  may  also  be  evaluated  by  forming  the  wedge 
product  of  the  row  vectors.  The  proof  of  these  theorems  will  be  found 
in  [2]. 


The  definition  of  the  determinant  by  means  of  a  wedge  product  is 

not  very  useful  for  calculating  the  value  of  a  determinant,  but  it  is 

useful  in  obtaining  theoretical  results  about  determinants.  For  example, 

we  shall  use  this  definition  to  obtain  Cramer's  rule  for  the  solution 

of  a  system  of  linear  equations.  Suppose  that  a-^a^, . . . ,an,b  are  given 

n-vectors  and  that  E.. £  are  n  unknowns  to  be  found  from  the 

12  n 

system  of  equations 


5lal  +  ?2a2  +' 


•+  5  a  =  b 
n  n 


(7) 


Let  us  solve  for  4^.  Since  the  wedge  product  of  a  vector  with  itself 
is  zero,  we  may  eliminate  by  multiplying  equation  (7)  with  a^. 

In  this  way  we  get 

*1  Va2  +  53  a3Aa2  +"*+  ?n  anAa2  =  bAV 
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Similarly,  multiplying  by  a^,...,an  will  eliminate  the  other  unknowns 
and  we  find 


or 


h  alAa2  A"‘A  an  "  bAa2  A‘“A  V 

bAa  a* • • a  a 

^  _  2 _ n_ 

1  a. Aa  a**’A  a  ’ 

12  n 


which  is  the  well-known  Cramer's  formula. 


Derivations 


Given  an  n-dimensional  linear  space  such  as  the  space  of  n-component 
vectors,  we  may  form  bi vectors  by  taking  the  wedge  product  of  two  vectors, 
tri vectors  by  taking  the  wedge  product  of  three  vectors,  and  so  on.  We 
shall  say  that  the  original  elements  of  the  space  are  elements  of  degree 
one,  bivectors  are  elements  of  degree  two,  tri vectors  are  of  degree  three, 
and  so  on  until  we  get  to  elements  of  degree  n.  We  may  form  the  sum  of 
elements  of  different  degrees  in  a  purely  formal  fashion  thus: 

3el  +  e2  ~  2elAe2  +  6elAe2Ae4’ 

Note  that  this  sum  cannot  be  simplified.  Since  we  can  also  always  form 
the  wedge  product  of  any  two  elements  of  arbitrary  degree,  we  have  an 
algebra  called  the  exterior  or  Grass man  algebra  over  the  linear  space. 

An  operator  D  which  maps  elements  of  the  exterior  algebra  into 
elements  of  the  algebra  is  said  to  be  a  linear  operator  if 


D(ax  +  Py)  =  aDx  +  pDy 
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for  any  elements  x  and  y  of  the  algebra  and  any  scalars  a  and  j3. 

If  the  operator  D  also  has  the  property  that  it  maps  elements  of  degree 
k  into  elements  of  degree  k  +  v  and  if  D  behaves  like  a  derivative 
with  respect  to  products,  that  is,  if 

D(xAy)  =  (Dx) Ay  +  (-)kV  xA(Dy)  (8) 

for  x  an  element  of  degree  k,  then  D  is  called  a  derivation  of  degree 
v.  For  example,  in  the  space  of  three-dimensional  vectors,  if  the  vectors 
are  functions  of  the  time  t,  then  ^  is  a  derivation  of  degree  zero. 
Another  example  is  given  by  the  operator  div  or  V*  .  This  operator 
is  a  derivation  of  degree  minus  one  because  it  maps  vectors  which  are  of 
degree  one  into  scalars  which  are  of  degree  zero.  The  product  rule  (8) 
is  also  obeyed  as  can  be  seen  from  the  well-known  formula 

V  •  (a  x  b)  =  (V  x  a) •  b  -  a  •  (V  *  b) . 

Notice  that  we  have  used  here  the  vector  product  instead  of  the  wedge 
product.  This  is  only  possible  in  three-dimensional  space. 

We  shall  prove  a  result  which  will  be  useful  to  us  in  the  next  chapter. 

2 

Theorem.  If  D  is  a  derivation  of  odd  degree  v,  then  D  is  a 
derivation  of  degree  2v . 

The  proof  is  obtained  by  straightforward  calculation.  We  have,  by 
definition, 


D(xAy)  =  (Dx) Ay  +  (-)kv  xA(Dy) 
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if  x  is  an  element  of  degree  k.  Then  since  Dx  is  an  element  of 
degree  k  +  v,  we  get 

D2(xAy)  =  (D2x) Ay  +  (-) v (Dx) A(Dy)  +  (-)kV[ (Dx) A(Dy)  +  (-)kv  xA(D2y)l 

=  (D2x) Ay  +  (-)2kV  xA(D2y), 

2 

because  (-)  +  (-)  =0  since  v  is  odd.  Thus  satisfies 

the  product  rule  (8)  and  since  it  maps  elements  of  degree  k  into 
elements  of  degree  k  +  2 v  it  is  a  derivation  of  degree  2v . 

Two  final  remarks.  First,  it  is  sufficient  to  define  a  derivation 
on  elements  of  degree  one  because  by  the  product  rule  (8)  it  can  be 
extended  successively  to  elements  of  arbitrary  degree.  Second,  no  matter 
what  the  degree  of  a  derivation  D  is,  if  D  is  applied  to  a  scalar,  it 
gives  zero.  This  follows  immediately  from  (8)  if  we  put  y  =  a,  a  scalar. 

We  find 

lev 

D(ax)  -  aDx  +  (-)  x  a  Da, 
or 

x  A  D(a)  =  0, 

for  all  a.  Therefore,  D(a)  =  0. 
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Chapter  II  Topology  of  Networks 


Introduction 

This  chapter  will  study  the  geometry  of  networks  with  no  reference 
to  their  electrical  properties .  We  shall  assume  that  the  network  is  a 
connected  set  containing  b  branches  and  n  +  1  nodes.  We  shall  assume 
that  each  branch  connects  two  nodes  and  that  each  node  lies  on  at  least 
two  branches  (see  Figure  l) .  Further  we  shall  assume  that  each  branch 


P„ 

i 

Figure  1 


carries  a  current  and  a  voltage  and  finally  we  shall  assume  that  these 
currents  and  voltages  satisfy  Kirchoff's  Laws,  namely,  the  sum  of  the 
currents  at  any  node  point  is  zero  and  the  sum  of  the  voltages  in  any 
closed  circuit  is  also  zero.  In  this  chapter  we  shall  assume  no  relation¬ 
ship  between  current  and  voltage.  These  relationships  will  be  considered 
in  the  next  chapter  where  we  will  discuss  the  electrical  properties  of 
networks . 

Chains 

Consider  a  network  of  b  branches  and  of  n  +  1  nodes.  In  Figure  1, 
we  have  an  example  of  such  a  network  with  b  -  8,  and  n  =  4.  We  shall 
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name  the  nodes  P„,P,  ,P„, . . .  ,P  .  The  nodes  will  be  called  O-cells.  A 
O’  1’  2  n  - 

linear  combination  of  O-cells  such  as 


5P, 


6P1  +  P2  "  P4 


will  be  called  a  0 -chain.  Later  we  shall  give  a  physical  interpretation 
for  O-chains.  At  present,  a  O-chain  will  be  just  a  formal  sum  of  the 
symbols  Pq>P]_» • •  •  >Pn *  The  set  of  all  O-chains  forms  a  linear  space. 


The  branch  joining  P.  and  P.  will  be  denoted  by  P.P.  if 

J  t  j 


the 


branch  is  to  be  considered  oriented  in  the  direction  from  P.  to  P . . 

i  J 

Since  the  direction  from  P^.  to  P^  is  the  opposite  of  the  direction 

from  P.  to  P.,  we  shall  write 
1  J 


P.P.  =  -  P.P. . 

1  J  J  i 


(1) 


In  this  way,  an  oriented  branch  may  be  considered  an  anti -commutative 
product  of  nodes . 


A  branch  will  be  called  a  1-cell  and  a  linear  combination  of  1-cells 
such  as 


2P  P 
0  1 


6  P2P1  + 


5P4P3 


will  be  called  a  1 -chain.  The  set  of  all  1-chains  forms  a  linear  space. 
We  shall  assume  that  a  chain  contains  only  those  1-cells  which  are 
branches  of  the  given  network.  Thus  in  Figure  1  no  chain  would  contain 
a  1-cell  such  as  PqP3  ‘  Tf  we  wan^  emPhasize  that  a  1-chain  contains 
only  branches  of  the  network,  we  shall  call  it  an  admissable  chain. 

A  closed  circuit  on  a  network  will  be  called  a  2 -cel  1 .  Thus,  in 
Figure  1,  P^P^Pg,  T?^T^Pq,  P^^P^  are  examples  of  2-cells.  Each  2-cell 
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will  be  considered  orientated  according  to  the  order  in  which  the  points 
are  described.  Thus  is  oriented  clockwise  and  Pq^2^1  is  oriented 

counter-clockwise.  We  shall  write 


PPP  =  -  P  P  P  . 
0  1  2  0  2  1 


Notice  that  this  equation  is  just  what  would  be  expected  if  the  2-cells 
are  considered  as  triple  products  of  nodes  with  the  products  obeying  the 
anti-commutative  law  (l). 


A  linear  combination  of  2-cells  such  as  -  P^P-Pq  -  3 PqP^P^ 

will  be  called  a  2-chain.  It  is  convenient  to  restrict  the  concept  of 

2-cells  to  symbols  containing  three  points  only;  consequently  the  circuit 

P-.P^P^P.  in  Figure  1  will  not  be  considered  as  a  2-cell  but  as  the  2-chain 
1  <■  j  U 

P-.P„P-  +  P-,P„P;.  Because  of  the  orientation  of  the  2-cells,  the  sum  in 
1  i  1  3  A 

the  2-chain  contains  the  common  side  P-^P^  traversed  twice,  once  in  the 

positive  and  once  in  the  negative  direction;  consequently,  the  chain  reduces 

to  the  branches  F,P„,  P-P.,,  P,P,,  P,Pn  traversed  in  that  order,  thus 
1  j  3  4  4-1- 

giving  the  circuit  P^PgP^P^. 

Boundary 

We  shall  define  an  algebraic  operator  on  the  spaces  of  chains  of  all 
dimensions.  This  operator  will  be  called  the  boundary  operator,  because 
when  this  operator  is  applied  to  cells  it  will  produce  the  geometrical 
boundary  of  the  cell.  We  denote  this  operator  by  the  symbol  3  and 
we  define  its  action  as  follows : 


The  operator  3  will  be  a  derivation  on  the  space  of  all  chains. 


When  it  is  applied  to  a  scalar,  it  gives  zero  as  any  derivation  will 
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When  3  is  applied  to  a  node,  it  gives  one.  If  we  now  consider  O-cells 
as  elements  of  degree  one,  1-cells  as  elements  of  degree  two,  and  2-cells 
as  elements  of  degree  three,  we  may,  by  the  remark  at  the  end  of  Chapter  I, 
extend  3  as  a  derivation  from  O-cells  first  to  1-cells  and  then  to  2-cells. 
In  doing  so,  we  shall  use  equation  (8)  of  the  preceding  chapter  and  we 
shall  assume  that  3  is  a  derivation  of  degree  minus  one. 


Let  us  consider 


3(P.  P.)  . 


We  have 


3(P.P.)  =  3(P.)P.  -  P.3(P.)  =  P.  -  P.  .  (2) 

i  J  i  J  -l  J  J  i 

We  see  that  the  boundary  of  the  1-cell  P^P^.  is  the  0 -chain  P^  -  P^, 

that  is ,  the  end  point  P .  counted  positively  and  the  initial  point  P. 

J  1 

counted  negatively.  Note  that  the  boundary  of  P_.P^  is  P^  -  P^ ,  which 
is  the  negative  of  the  boundary  of  P^P^. . 


Since  3  is  a  linear  operator,  the  boundary  of  a  1-chain  may  be 
found  by  applying  3  to  each  term  in  the  chain.  Thus 

3(PiP2  +  2P0P1)  2Pj.  -  2P0  =  P2  +  Pj  -  2P0- 


Let  us  now  consider  the  boundary  of  the  2-cell  P.P.P,  .  Again  using 

1  J  k 

(8)  of  the  preceding  chapter  and  (l)  and  (2)  of  this  chapter,  we  find 
that 


a(p.p.p.)  =  a(p.)p.p,  -  p.  3(p .p,  )  =  p.p,  -  p.p,  +  p.p.. 

v  l  j  k'  '  l'  J  k  l  v  J  k'  j  k  xk  ij 


We  see  that  the  boundary  of  the  2-cell  P.P.P,  in  Figure  2  is  the  1-chain 

1  .1  & 


.P.  +  P.-F.  -  P.P- 
l  3  J  k  l  . 
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that  is,  the  branches  that  form  the  boundary  of  the  2-cell  P^P.P^,  each 


P. 

1 


Figure  2 


counted  with  a  plus  or  minus  sign  according  as  the  orientation  of  the 
branch  agrees  or  does  not  agree  with  the  orientation  of  the  2-cell. 


Again,  since  3  is  a  linear  operator,  the  boundary  of  an  arbitrary 
2-chain  may  be  found  by  applying  3  to  each  term  in  the  chain.  Thus  the 
boundary  of  P-^P^P^  +  P-^P^P^  is 


P1P2  +  P2P3  '  P1P3  +  P1P3  +  P3P4~  P1P4 


P1P2  +  P2P3  +  P3P4  "  P1P4* 


Notice  that  this  1-chain  actually  represents  the  geometrical  boundary 

of  the  circuit  P,P  P„P,  in  Figure  1. 

X  2  3  4 


We  shall  now  prove  a  fundamental  topological  fact,  namely, 


Theorem.  The  boundary  of  a  boundary  is  zero. 


For  the  proof,  we  use  the  theorem  we  proved  at  the  end  of  the 

2 

preceding  chapter.  Since  9  is  a  derivation  of  odd  degree,  3  is  a 

2 

derivation  of  degree  minus  two,  but  3  (P. )  =  3(l)  =  0.  Therefore, 

2 

when  we  use  the  remark  at  the  end  of  Chapter  I,  we  conclude  that  3 
extended  to  all  chains  gives  zero. 
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Current  Chains 


A  chain  may  be  interpreted  as  currents  carried  in  the  appropriate 


cells.  For  example,  a  0-chain  such  as  Sa^P^.  will  describe  currents 

of  magnitude  flowing  into  the  node  P^.  We  shall  call  such  a  chain 

a  node-current  chain.  A  1-chain  such  as  ZB.  .P.P.  will  describe  currents 
-  l  3 

of  magnitude  flowing  along  the  branch  P^P^..  We  call  such  a  chain 

a  branch-current  chain.  A  2-chain  such  as  Zy.  .nP.P.Pn  will  describe 
currents  of  magnitude  flowing  in  the  circuit  P^P^P^.  We  call 

such  a  chain  a  loop -current  chain. 


Consider  a  branch  current  aP.P,  .  Its  boundary  is  aP  -  aP.  which 

3  k  k  j 

may  be  interpreted  as  a  current  of  magnitude  a  going  into  the  node  P^ 

and  a  current  of  magnitude  o  going  out  of  the  node  P..  Generalizing, 

3 

we  see  that  the  boundary  of  a  branch-current  chain  is  the  node-current 
chain  which  describes  the  total  amount  of  current  flowing  into  the  nodes 
because  of  the  branch  currents . 


We  can  obtain  a  similar  result  for  loop-current  chains .  The 
boundary  of  a  loop-current  chain  is  the  branch-current  chain  which 
describes  the  currents  flowing  through  the  branches  that  form  the  sides 
of  the  loop. 

Suppose  we  have  a  branch -current  chain  which  is  the  boundary  of  a 
loop-current  chain.  Using  the  theorem  about  the  boundary  of  a  boundary 
being  zero,  we  find  that  the  boundary  of  this  branch -current  chain  must 
be  a  node-current  chain  which  is  identically  zero.  This  means  that  the 
sum  of  the  currents  flowing  into  any  node  must  be  zero;  consequently, 
Kirchoff's  Node  Law  is  satisfied.  Conversely,  if  the  network  is  connected, 
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Kirchoff's  Node  Law  will  imply  that  the  branch -current  chain  is  the 
boundary  of  a  loop-current  chain.  We  shall  state  this  result  in  the 
following  way :  A  branch -current  chain  must  be  a  boundary  1-chain. 

Co -chains 

To  discuss  voltages  in  the  same  way  we  have  discussed  currents,  we 

introduce  the  concept  of  a  co-chain  and  a  co-boundary.  A  k-co-chain 

is  a  numerical -valued  function  defined  on  k-cells.  For  example,  consider 

the  O-co-chain  c^  defined  as  follows:  The  value  of  c^  for  any 

P.  j  /  2,  is  zero;  the  value  of  c  for  P„  is  one.  Note  that  this 
J  !  U  *- 

co— chain  may  be  extended  linearly  to  be  a  function  of  chains;  thus, 

C0(2akPk)  =  V 


(4) 
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For  example,  the  co-boundary  of  the  O-co-chain  defined  previously  is 

the  1-co-chain  whose  value  on  the  1-cell  P. .  is  zero  if  neither  i 

iJ 

nor  j  is  two;  on  the  1-cell  P.  _  its  value  is  1;  and  on  the  1-cell 

1,2 

P„  .  its  value  is  -  1.  Notice  that  the  co-boundary  is  zero  if  P 
2,  J  2 

is  not  on  the. boundary  of  the  1-cell,  but  it  is  +1  if  is  the  end 

point  of  the  1-cell  and  -  1  if  it  is  the  origin  of  the  1-cell. 

We  have  proved  that  the  boundary  of  a  boundary  is  zero.  Using  this 
result,  we  can  prove  the  following  dual  theorem: 


The  co-boundary  of  a  co-boundary  is  zero. 


The  proof  is  a  simple  consequence  of  (4)  •  We  have,  for  any  p-chain 


dP, 


,2 

6  c 


P-2 


dp  -  5c 


P-2 


9dP  =  c 


P-2 


92dp  =  0 


2  2 

because  of  the  fact  that  5  =0.  Since  the  value  of  6  c  2  on  any 

2 

p-cham  is  zero,  we  must  have  5  c  _  =  0. 

r  ’  p-2 


Voltage  co-chains 


Just  as  currents  were  defined  by  chains,  so  voltages  will  be  defined 
by  co-chains .  A  0-co-chain  will  define  the  value  of  a  potential  at 
each  node;  a  1-co-chain  will  define  a  voltage  in  each  branch;  a  2-co¬ 
chain  will  define  a  total  voltage  for  a  2-cell.  For  example,  consider 
the  0-co-chain  which  states  that  the  potential  at  P^  is  2,  at  P^ 
is  4,  and  at  P^  is  5«  We  shall  write  this  as 

2pQ  +  4?^^  +  5p2- 
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Note  that  for  k  =  0,1,2  is  actually  the  co-chain  that  has  the 

value  1  for  and  the  value  zero  for  any  other  0-cell.  Let  us 

consider  the  co-boundary  of  this  co-chain.  It  is 

2poPi  +  3PqP2  +  P1P2* 

Where  p.p.  is  the  co-chain  that  has  the  value  1  for  the  1-cell 
i  J 

P.P.,  the  value  -  1  for  the  1-cell  P.P.  and  the  value  zero  for 
i  J  3  i 

any  other  1-cell.  We  see  that  the  co-boundary  of  a  node-potential 
co-chain  is  the  branch -voltage  co-chain  such  that  the  voltage  in  each 
branch  is  the  difference  of  the  node  potentials  at  the  ends  of  the  branch. 

A  similar  calculation  will  show  that  the  co-boundary  of  a  branch- 
voltage  co-chain  is  a  loop-voltage  co-chain  such  that  the  voltage  in  each 
loop  is  the  oriented  stun  of  the  voltages  in  the  branches  bounding  the 
loop. 


Suppose  we  have  a  branch-voltage  co-chain  which  is  the  co-boundary 
of  a  node-potential  co-chain.  Since  the  co-boundary  of  a  co-boundary 
is  zero,  we  see  that  the  co-boundary  of  the  branch-voltage  co-chain 
must  be  zero.  This  result  implies  that  the  sum  of  the  voltages  in  any 
loop  must  be  zero.  Thus,  Kirchoff's  Circuit  Law  is  satisfied.  Conversely, 
it  can  be  shown  that  if  the  network  is  connected  and  if  Kirchoff's  Law 
is  satisfied,  then  the  branch -voltage  co-chain  must  be  the  co-bcsundary 
of  a  node -potential  co-chain.  We  shall  state  this  result  as  follows: 

The  branch  voltages  must  form  a  co-boundary  co-chain. 

Tellegen's  Theorem 

We  shall  prove  an  important  theorem  connecting  the  branch  currents 
and  the  branch  voltages  in  an  arbitrary  network  if  Kirchoff's  Laws  are 
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satisfied.  The  theorem,  first  proved  by  Tellegen,  is  as  follows: 

Let  I  be  a  current  vector,  that  is,  let  I  be  a  b-component 
vector  such  that  each  component  represents  the  current  in  a  particular 
branch  of  the  network.  Let  V  be  a  voltage  vector,  that  is,  let  V 
be  a  b-component  vector  such  that  each  component  represents  the  voltage 
in  a  branch.  If  Kirchoff's  Laws  are  satisfied,  then  the  scalar  product 
of  I  and  V  must  be  zero. 

Notice  that  no  relation  between  I  and  V  is  assumed.  All  that 
is  required  is  that  the  current  and  voltage  vector  each  separately  satisfy 
Kirchoff's  Laws.  The  proof  of  the  theorem  will  follow  easily  from  the 
topological  characterization  of  branch-currents  and  branch-voltages  given 
previously.  To  the  vector  V  there  corresponds  a  co-boundary  1-co-chain 
m  C-^  and  to  the  vector  I  there  corresponds  a  boundary  1-chain  . 

From  (3)  we  see  that  the  scalar  product  of  I  and  V  is  equal  to  c^  •  d^ . 
Since  c^  is  a  co-boundary,  it  may  be  written  as  6c ^  and  since  d^ 
is  a  boundary,  it  may  be  written  as  9d  .  Now  by  (4) . 

6cq  •  3d2  =  c Q  •  32d2  =  0, 

because  the  boundary  of  a  boundary  is  zero;  therefore  the  theorem  is 
proved . 

Tellegen's  theorem  can  be  generalized  and  written  in  an  integral 
form.  Suppose  we  consider  the  behavior  of  a  network  as  a  function  of 
some  parameter  such  as  the  time  t.  For  each  value  of  t  in  some  range 
r  there  will  be  a  current  vector  l(t)  and  a  voltage  vector  V(t) . 

We  shall  prove  the  following: 
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If  I  and  V  independently  satisfy  Kirchoff's  Laws,  then 
5  Vdl  =  0. 

r 

The  proof  is  an  easy  consequence  of  Tellegen's  theorem.  Since 
l(t  +  h)  and  I(t)  are  both  current  vectors,  we  have 

V(t) [I(t  +  h)  -  I(t)]  =  0 

for  all  values  of  t  and  h.  Letting  h  go  to  zero,  we  find  that 
V(t)dl(t)  =  0 

for  all  values  of  t.  Approximating  the  integral  in  the  theorem  by  a 
sum,  we  finally  obtain  the  result 

5  Vdl  =  0. 

r 

Notice  that  a  similar  argument  will  prove  the  dual  result,  namely. 


j  IdV  =  0. 

r 
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Chapter  III  Non-linear  Networks 
Relations  between  currents  and  voltages 

In  the  preceding  chapter  we  have  considered  the  branch  currents 
and  the  branch  voltages  as  completely  independent  entities,  limited  only 
by  Kirchoff’s  Laws.  However,  in  any  network  the  current  in  a  branch  is 
related  to  the  voltage  in  that  branch  by  an  electrical  device  such  as  a 
resistor,  an  inductor,  or  a  capacitator.  Kirchoff 's  Laws  together  with 
the  electrical  relations  between  the  branch  currents  and  the  branch 
voltages  will  give  us  a  complete  set  of  equations  from  which  the  elec¬ 
trical  behavior  of  the  network  can  be  determined. 

We  begin  by  specifying  the  relation  between  the  orientation  of  the 
branch  and  the  signs  of  the  branch  current  and  the  branch  voltage.  If 
the  current  through  a  branch  is  in  the  same  direction  as  the  orienta¬ 
tion,  the  current  will  be  called  positive;  otherwise,  the  current  will 
be  called  negative.  The  voltage  across  a  branch  will  be  equal  to  the 
potential  at  the  endpoint  of  the  branch  minus  the  potential  at  the 
initial  point  of  the  branch.  Since  the  current  goes  from  a  node  of 
higher  potential  to  one  of  lower  potential,  our  conventions  about  signs 
will  produce  an  extra  minus  sign  in  the  customary  relations  between 
current  and  voltage.  For  example,  consider  Figure  1  in  which  the  upper 
node  of  the  branch  is  at  six  volts  potential  above  the  potential  of  the 
lower  node  and  in  which  the  resistance  is  two  ohms.  Of  course,  the 
current  through  the  branch  is  three  amperes  flowing  downward.  If  the 
branch  is  oriented  downward,  then  the  current  is  positive  and  the 
voltage  is  negative.  However,  if  the  branch  is  oriented  upward,  then 
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Figure  1 

the  current  is  negative  and  the  voltage  is  positive.  In  both  cases, 
the  relation  between  the  voltage  v  across  the  branch  and  the  current 
i  through  it  is  v  =  -  Ri,  where  R  is  the  resistance  of  the  branch. 

In  addition  to  linear  resistors,  we  shall  also  consider  non-linear 
resistors  in  which  the  relation  between  current  and  voltage  is  given 
by  the  equation  f(i,v)  =  0.  This  equation  can  be  solved  for  either  of 
the  variables  i  or  v  as  a  function  of  the  other  variable.  We  shall 
not  require  that  either  of  the  resulting  functions  be  single-valued 
because  we  shall  permit  the  characteristic  of  the  resistor  to  have  the 
shape  shown  in  Figure  2. 


Figure  2 
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The  inductors  and  capacitors  may  also  be  non-linear.  For  an 
inductor,  the  relation  between  current  and  voltage  may  have  the  form 

•v 

where  L(i)  is  a  non-negative  inductance.  For  a  capacitator,  the 
relation  between  current  and  voltage  may  have  the  form 

i  =  -  C( v)  g 

where  C(v)  is  a  non -negative  capacitance. 

Equilibrium  equations 

Let  us  consider  again  a  connected  network  containing  n  +  1  nodes 
and  b  branches.  A  set  of  branches  which  is  connected  and  contains  no 
closed  loop  is  called  a  tree .  A  maximal  tree  is  a  tree  which  is  not 
contained  in  any  other  tree.  It  is  easy  to  see  that  a  maximal  tree  must 
contain  all  the  nodes;  for,  if  a  node  did  not  belong  to  the  tree,  a 
branch  joining  this  node  to  the  tree  could  be  added  without  spoiling  the 
tree  shape.  It  is  also  easy  to  see  that  a  maximal  tree  must  have  exactly 

n  branches;  for,  starting  with  one  fixed  node,  we  must  add  one  branch 

for  each  node  that  is  to  be  connected  to  this  fixed  node. 

The  network  will  have  many  maximal  trees .  Choose  one  such  tree  and 
consider  the  b  -  n  =  L  branches  that  do  not  belong  to  the  tree.  These 

L  branches  are  called  links  of  the  tree.  Note  that  if  any  link  is  added 

to  a  maximal  tree  the  resulting  network  will  have  a  loop  in  it. 

From  the  definition  of  a  tree  it  is  clear  physically  that  the 
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voltages  in  the  n  branches  of  a  maximal  tree  can  be  chosen  arbitrarily 
because  the  tree  has  no  loops.  Then  by  adding  each  link  in  turn  to  the 
tree  and  using  Kirchoff 's  Circuit  Law  we  may  find  the  voltage  in  that 
link.  Since  the  network  contains  L  links,  this  argument  shows  that 
there  are  L  independent  loop  equations . 

Instead  of  starting  with  voltages,  we  may  begin  with  currents.  As 
we  add  a  link  to  the  maximal  tree,  we  may  impress  an  arbitrary  current 
on  it.  It  is  clear  that  when  all  the  links  with  their  currents  have 
been  added  to  the  tree ,  the  currents  in  the  other  n  branches  of  the 
tree  can  be  determined.  This  argument  implies  that  there  are  exactly 
n  independent  mesh  equations .  Notice  that  so  far  we  have  not  used  the 
electrical  properties  of  the  branches  but  have  used  only  Kirchoff's 
Laws . 


Suppose  we  pick  a  set  of  r  <  L  branches  whose  currents  we  shall 

call  i^,...,ir  and  a  set  of  s  <  n  branches  whose  voltages  we  shall 

call  v  ....... v  ,  .  We  assume  that  no  branch  has  both  a  current  and 

r+1’  r+s 

a  voltage  specified.  We  shall  show  that  all  the  other  currents  and 
voltages  can  be  determined  from  these  specified  currents  and  voltages. 

First,  from  Kirchoff's  Laws  we  have  L  equations  for  the  voltages  and 
n  equations  for  the  currents.  Then,  in  the  b-r-s  branches  for  which 
neither  current  nor  voltage  was  specified,  we  use  the  electrical  properties 
of  these  branches  to  obtain  b-r-s  equations.  Thus,  we  have  a  total  of 
L+n+b-r-s=  2b-r-s  equations  to  find  the  unspecified  b-r  currents 
end  the  b-s  voltages .  Since  the  number  of  equations  is  equal  to  the 
number  of  unspecified  quantities,  we  assume  that  they  can  be  solved  to 
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give  equations  of  the  form 


i 


v 

P 


=  G  (i,v)  • 

r 


1  <  (i  <  b. 


(1) 


Notice  that  in  this  discussion  we  have  not  used  the  electrical 
properties  of  the  branches  in  which  either  the  current  or  the  voltage 
was  specified.  If  we  assume  that  the  voltage  in  the  r  links  is 
related  to  the  current  by  equations  of  the  form 

vn  =  >  1  <  P  <  r 

P  P  P 

and  if  we  assume  that  the  current  in  the  s  branches  is  related  to  the 
voltage  by  equations  of  the  form 


i  =  f  (v  )  ,  r  +  1  <  <3  <  r  +  s 

d  d  d  — 


we  obtain  the  following  r  +  s  equations  from  which  the  equilibrium 
state  of  the  network  can  be  determined: 


g  (i  )  =  Gp(i»v) ,  1  <  P  <  r, 

f  (v  )  ~  F  (i,v)  >  r+l<c<r+s. 

do  d  —  — 


(2) 


An  example  may  clarify  the  argument.  Consider  the  network  in  Figure  3 
for  which  n  =  2  and  b  =  A.  Branches  numbered  3  and  4  will  form  a  maxi¬ 
mal  tree  and  branches  1  and  2  will  be  links  for  this  tree.  Let  us  assume 
that  the  currents,  i„  and  i^  in  branches  1  and  2  are  specified  and 
also  the  voltage  v^  in  branch  3*  This  implies  that 
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r  =  2  and  s  =  1.  By  the  use  of  Kirchoff 's  Laws  we  find  that  the 
currents  in  branches  3  and  4  are  given  by  the  equations  i^  =  i^  +  i^, 
i.  =  i-  +  i„.  Next,  assuming  that  the  electrical  properties  of  the 
unspecified  branch  4  are  such  that  the  relation  between  current  and 
voltage  is  given  by  v  =  g.(i)  we  find  that 


and 


v4  =  g4(i4}  =  g4(Il  +  I2) 


V1  V2 


V4  "  v3  =  "  v3  "  g4<‘il  +  i2)  ‘ 


Thus,  all  the  unspecified  variables  have  been  expressed  in  terms  of  the 
specified  ones.  Now,  using  the  electrical  properties  of  toe  specified 
branches,  we  finally  obtain  the  following  equilibrium  equations  for  the 
network: 


gl^l^  s2<'12>  ~  V3  g4^1l  +  12> 


f3(v3)  =  +  i1  +  i2. 


In  writing  this  set  of  equilibrium  equations  we  have  assumed  that  all 
the  branches  were  resistors.  Suppose,  however,  that  branches  1  and  2 
were  inductances  and  that  branch  3  was  a  capacitor.  Then  the  dynamical 
equations  for  the  network  would  be 


dt  L2^12^  dt  +  V3  +  g4^Xl  +  12^  ’ 
_  dv  „ 

”  C3(v3}  dF  _  X1  +  x2' 


(3) 


Notice  that  any  set  of  currents  and  voltages  can  be  specified,  as 
long  as  Kirchoff 's  Laws  are  not  violated.  The  only  reason  for  specifying 
a  particular  set  of  currents  and  voltages  is  that  we  wish  to  write  the 
equilibrium  equations  of  the  network  in  terms  of  this  particular  set. 
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Consider  the  two  terms  in  the  bracket  on  the  right-hand  side.  Using  (4) 

for  r  +  1  ^  p  <  s,  we  can  express  the  sum  in  the  bracket  as  a  function 

of  the  specified  variables  i,v  at  the  end  of  r  and  at  the  beginning 

of  P.  For  the  integral  term  in  the  bracket,  notice  that  the  branches 

involved  are  only  the  resistive  ones  in  which  v^  depends  only  on  i^. 

(it  is  immaterial  that  the  relation  between  i  and  v  may  not  be 

single-valued.  The  initial  point  of  r  determines  a  point  on  the  curve 

f  (i  v  )  =  0  and  we  just  integrate  along  that  curve.)  The  point  of  the 
P  P  P 

argument  is  that  a  term  such  as 


f  v  di  =  f  g  (i  )di 

J  1 1  J  \  ii' 


P  p 


P  P 


depends  only  on  the  endpoint  of  P  and  by  (4)  this  is  a  function  of  the 
specified  variables  only.  Our  conclusion  is  that  the  bracket  term  in 
(5)  is  a  function  P(i,v)  of  the  specified  variables  only.  We  call  this 
function  the  mixed  potential  of  the  network. 


We  may  then  write  (5)  as  follows: 


r  r+s 

2  s  vpdip  -  2  s  ipdvp  +  p^)  =  °-  (6) 

i  r  h  *  r+i  r  ^  p 


Since  the  variables  i,v  are  independent,  we  see  that  (6)  implies  that 


9P 


3P 

3v 

P 


1  <  p  <  r 

r  +  l<p<r  +  s. 


But  v  ,  1  <  p 
r  +  1  <  p  <  s 


<  r,  is  the  voltage  across  an  inductor  and 
is  the  current  across  a  capacitator.  Using 


the 
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electrical  relations  for  these  branches,  we  finally  obtain  the  dynamical 
equations  of  the  network,  namely 


di 

t  ( ■  \  ^  -  dp 

p^p  dt  di  ’ 

p, 


r  (  ap 

W  at  -  dV 


1  <  P  <  r 


r  +  l<p<r  +  s, 


(7) 


where 


P(i.v) 


b 

X 


5  v,,di,,  + 


r+s+1  r 


P  p 


r+s 

Z 

r+1 


i  v 


(8) 


For  convenience,  we  have  dropped  the  tilde  signs. 

The  Mixed  Potential  for  a  Complete  Set 

The  mixed  potential  can  be  given  a  particularly  useful  form  in  case 
the  r  currents  through  the  inductors  and  the  s  voltages  through  the 
capacitators  form  a  complete  set  of  variables.  This  means  they  are  a 
set  of  currents  and  voltages  which  do  not  violate  Kirchoff's  Laws  and 
which  are  such  that  either  the  current  or  the  voltage,  or  both,  in  every 
branch  is  determined  by  Kirchoff's  Laws  only.  For  example,  and 

v^  form  a  complete  set  of  variables  for  the  network  in  Figure  3  because, 
no  matter  what  values  they  have,  they  do  not  violate  Kirchoff's  Laws 
because  the  current  in  the  unspecified  branch  4  is  determined  by  Kirchoff's 
Laws  to  be  i-^  +  i 

Consider  a  network  in  which  the  specified  variables  form  a  complete 
set.  Let  Nv  be  the  set  of  branches  of  the  network  in  which  the  voltages 
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are  specified  by  Kirchoff's  Laws  and  let  fL  be  the  complementary  set 
of  branches,  that  is,  those  in  which  the  current  is  specified  by  Kirchoff's 
Laws .  For  example ,  in  Figure  3 ,  N  would  contain  only  the  branch  3  and 
IL  would  contain  branches  1,2,  and  4.  Using  this  decomposition  of  the 
network,  we  may  write  (8)  as  follows: 


P(i.v) 


r+s 

=  [  2  +  Z  ]  J  v  di  +  Z  i  v  I  . 

H>r+s  n>r+s  r  ^  ^  r+1  ^  11  'r 

peNL  ^eNv 


The  second  sum  in  the  bracket  can  be  integrated  by  parts  to  give 


P(i,v)  =  2  J 

p>r+s  r 


v  di 

p,  p 


p>r+s 


i  dv 


2  iv 
peNy  ^  ^'r 


(9) 


Notice  we  have  used  the  fact  that  the  inductive  branches  belong  to  N. 
and  the  capacitative  branches  belong  to  Nv  to  simplify  one  of  the  sums 
in  (9). 


The  individual  terms  in  (9)  can  be  given  an  interesting  physical 
interpretation.  The  first  term 

5  v  di 
f  P  P- 

is  well-defined  as  a  line  integral  over  the  curve  connecting  i  and  v 
in  the  branch  p.  We  call  this  integral  the  current  potential  of  the 
branch  p.  Similarly,  the  integral 

5  i  dv 
P  P  P 

is  called  the  voltage  potential  of  the  branch  p.  Note  that  in  (9)  the 
current  potential  is  evaluated  only  for  branches  in  IL  and  the  voltage 
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potential  only  for  branches  in  Nv»  Since  the  current  is  specified  in 
1L  and  the  voltage  in  N  ,  we  see  that 

P(i,v)  =  F(i)  -  G(v)  +  2  i  v  I  ,  (10) 

peNv  ^  ^  !r 

where  F(i)  is  the  current  potential  and  G(v)  is  the  voltage  potential 
It  can  be  shown  that  the  third  term  on  the  right-hand  side  of  (10)  is 
a  bilinear  function  of  the  specified  currents  and  voltages.  For  details 
we  refer  to  [l] . 

Example  of  the  Calculation  of  the  Mixed  Potential 

Consider  again  the  network  illustrated  in  Figure  3,  but  now  we  shall 
specify  its  electrical  properties  as  illustrated  in  Figure  4.  The  symbol 


in  branch  3  indicates  the  presence  of  a  non-linear  resistor  whose 
characteristic  will  be  assumed  to  be  of  the  form  v^  =  -  f (i^) . 

Since  there  is  only  one  inductor  and  one  capacitor  in  the  network 
of  Figure  4,  the  appropriate  variables  to  specify  the  dynamical  behavior 
of  this  system  are  i„,  the  current  across  the  inductor,  and  v.,  the 
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voltage  across  the  capacitor.  Note  that  i0  and  v.  do  not  form  a 

2  4 

complete  set  of  variables  because  the  use  of  Kirchoff's  Laws  alone  without 
the  electrical  properties  of  the  network  does  not  determine  the  value 
of  either  i?  or  v^. 


Nevertheless,  even  though  the  set  of  variables  is  not  complete,  a 
mixed  potential  does  exist  and  the  dynamical  equations  for  the  network 
in  Figure  4  can  be  written  in  the  form  (7) •  To  see  this,  let  us  write 
Kirchoff's  Laws  using  the  electrical  properties  of  the  network.  It  is 
clear  that  i«  =  -  i^  -  i^.  We  have 


di. 


L2  "dt"  = 


v4  =  f(ll  +  i2)  +  ilPd 


dv 


C4  dt"  "  W  +  12>  ’ 


(11) 


Since  we  have  three  equations  for  the  three  variables  i,  ,i  ,v, ,  the 

12  4 

system  is  completely  determined. 

Notice  that  the  second  equation  of  the  set  (ll)  is  not  a  differential 

equation.  This  occurrence  is  a  result  of  the  fact  that  the  set  of 

variables  i„  and  v.  is  not  a  complete  set.  We  can  of  course  eliminate 
2  4 

i^  from  the  equations  (ll)  and  thus  obtain  a  set  of  differential  equations 
only.  To  do  this,  we  notice  that  the  second  equation  of  (ll)  may  be 
written 


v4  +  12R1  =  f(il  +  i2)  +  (il  +  i2)Rl> 


(12) 
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so  that  v^  +  i^R^  ds  a  f^inc',:'ion  °f  +  i2*  ^et  us  ^nvel"t  this 
functional  relationship  to  get  i^  +  i^  as  a  function  of  v^  +  i^R^ 
and  let  us  write  this  new  relation  as 


11  +  i2  =  g(v4  +  i2Rl) 


New  eliminating  i^  from  the  first  and  third  equations  of  (ll),  we 
get  the  desired  result: 

dv4  _ 


C4  df  =  -  g(v4  +  i2Rl) 


dl2 

L2  "dt”  "  Rlg^v4  +  12R1^  "  ^2 


(13) 


The  set  of  equations  (13)  may  be  written  in  the  form  (7)  if  we 
introduce  the  mixed  potential  P  defined  as  follows: 

(V4+i2Rl) 

P  =  5  g( v)dv  -  jH.i'T.  (14) 

0 

Notice  that  the  form  (14)  for  the  mixed  potential  is  not  in  the  form  (10) , 
that  is ,  it  is  not  the  sum  of  a  term  depending  only  on  the  voltage  plus 
a  term  depending  only  on  the  current  plus  a  term  which  is  bilinear  in 
the  current  and  the  voltage.  Instead,  the  first  term  on  the  right-hand 
side  of  (14)  is  a  complicated  function  of  both  the  current  and  the  voltage. 

If  we  modify  the  network  in  Figure  4  so  that  it  is  the  network 
illustrated  in  Figure  5,  we  shall  have  an  example  of  a  network  containing 
a  complete  set  of  variables  and  we  shall  find  a  mixed  potential  of  the 
form  (10) .  Since  there  are  two  inductors  and  one  capacitor  in  Figure  5, 
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Equation  (16)  shows  that  P  has  the  form  (10)  with  the  current 


potential  equal  to 


1  R  ±2 
~  2  R1  1 


(il+i2) 

-  S  f(i)di, 

0 


no  voltage  potential  term,  and  the  remaining  part  being  the  bilinear 
function  v^(i^  +  i-2)  *  The  f°rmula  (16)  can  be  written  down  immediately 
from  (10)  without  going  through  (15) .  Prom  (9)  and  (10) ,  F(i)  is 
the  current  potential  evaluated  along  the  resistors  in  the  branches  in 
which  the  current  is  defined  by  Kirchoff's  Laws  alone  without  the 
electrical  properties  of  the  network.  In  Figure  5,  the  current  is 
defined  in  all  the  branches  but  since  there  are  resistors  only  in 
branches  1  and  3,  we  find 


F(i)  =  -  |  RjiJ 


-  5 

0 


(il+i2) 


f  (i)di. 


(17) 


The  voltage  potential  G(v)  must  be  evaluated  along  the  resistors  in 
the  branches  in  which  the  voltage  is  defined  by  Kirchoff's  Laws  alone. 
The  only  such  branch  is  4  and  since  this  branch  has  no  resistance,  we 
have  G(v)  =0. 


The  remaining  part  of  the  mixed  potential  is  defined  in  (9)  as 

Z  iaval  ’  U8) 

peN  ^ 
v 

where  the  sum  is  to  be  taken  over  all  branches  in  which  Kirchoff's  Laws 

specify  the  voltage.  Again,  the  only  such  branch  is  4>  and  since  the 

current  in  this  branch  is  -  (in  +  i„)  and  the  voltage  is  -  v.  because 

I  «c  4 

of  the  assigned  orientation  of  the  branch,  we  find  that  (18)  equals 
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v^(i^  +  i^) •  Combining  (17)  and  (18) ,  we  get 

(il+i2} 

P  =  -  2  Vl  -  *  f(i)di  +  V4(1l  +  ±2)> 

in  agreement  with  (16)  . 

Vector  Form  for  the  Dynamical  Equations 


Consider  a  network  with  a  complete  set  of  variables  i,,...,i  , 
Vr+1’ * ' * ’ Vr+s *  denote  this  set  of  variables  by  the  (r  +  s) -component 
vector  x.  Let  L  denote  the  diagonal  matrix  whose  entries  are  the 
inductances  L^(i) , . . . ,Lr(i)  and  let  C  denote  the  diagonal  matrix 
whose  entries  are  Cr+^(v)  > • • • >Cr+s(v) •  Pat 


where  the  zeros  denote  rectangular  matrices  whose  entries  are  all  zero. 
If  a  dot  denotes  differentiation  with  respect  to  time,  we  may  write  (7) 
in  the  following  vector  form: 

(19) 


-  J  x  = 


i-  =  apU) 

3x 


ap 


where  denotes  the  gradient  of  the  function  P. 


Consider  any  solution  of  the  equations  (19) .  This  solution  is  defined 
by  specifying  x  as  a  function  of  t;  say,  x  =  X(t).  Note  this  speci¬ 
fication  defines  a  curve  r.  Let  us  consider  the  way  in  which  the  mixed 
potential  varies  along  T.  We  have 

dP(X(  t))  _  3P  dx  _  /  ty 

a  +  Ay  ’  a+  <JX,X>  . 


(20) 


The  expression  on  the  right-hand  side  of  (20)  is  a  quadratic  form 
in  the  vector  X(t) .  Suppose  that  J  is  a  positive-definite  matrix, 
then  (20)  shows  that  the  value  of  the  mixed  potential  always  decreases 
along  any  curve  T.  Considering  the  mixed  potential  as  some  kind  of 
"energy"  of  the  network,  we  are  led  to  believe  that,  because  the  "energy" 
always  decreases,  the  network  must  always  be  stable.  This  semi -intuitive 
reasoning  will  be  made  precise  in  the  next  section  where  we  discuss 
Liapunov's  Stability  Theorem. 

Liapunov  Stability 

Let  x  denote  an  n-component  vector  and  let  F(x)  be  a  continuous 
and  differentiable  vector-valued  function  of  x.  Consider  the  vector 
system  of  first-order  equations 

X  =  F(X)  •  (21) 

The  equilibrium  points  of  this  system  are  the  points  x  in  n-space  E^ 
at  which  F(x)  =  0.  An  invariant  set  of  (2l)  is  a  set  of  points  in 
such  that  if  a  solution  of  (21)  starts  in  the  invariant  set  at  t  =  0, 
then  for  all  future  time  the  solution  remains  in  the  invariant  set.  It 
is  clear  that  the  equilibrium  points  of  (2l)  form  an  invariant  set. 

Let  V(xT  be  a  scalar-valued  function  of  x  with  continuous  first 
partial  derivatives.  Along  any  solution  of  (21)  specified  by  x  =  X(t), 
we  have 

*x(t»  •  grad  V  . 
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Let  us  put 

V(x)  =  F(x)  •  grad  V(x)  , 

so  that 

=  V(X(t)).  (22) 

If  V(x)  <  0  and  V(x)  >  0  for  some  region  of  containing 

the  origin,  then  V(x)  will  be  called  a  Liapunov  function.  We  shall 
prove  the  following  [3] : 

Theorem  . .  Let  R„  denote  the  region  around  the  origin  in  E.. 

u  n 

where  V(x)  <  C.  Suppose  that  R_  is  bounded  and  that  inside  R_ 

V(x)  >  0  for  x  ^  0,  and  V(x)  <  0.  Let  M  be  the  largest  invariant 
set  of  ( 21)  contained  in  the  set  of  points  for  which  V( x)  =  0 ;  then, 
every  solution  of  ( 21)  which  starts  in  R  at  t  =  0  tends  to  M  as_ 
t  goes  to  infinity. 

To  prove  the  theorem,  we  consider  any  solution  x(t)  of  (21)  which 

is  in  R  at  t  =  0;  therefore  V(x(0))  <  C.  Because  of  (22)  and  the 
u 

conditions  on  V(x)  ,  the  value  of  V  along  this  solution  is  non -increasing 
for  all  t  >  0,  and  V(X( t) )  <  C  for  all  t  >  0;  we  conclude  that  the 
solution  X(t)  can  never  reach  the  boundary  of  R  at  which  V(x)  =  C 
and  therefore  X(  t)  must  remain  in  R  .  Since  V(X(t))  is  non-increasing 

L» 

and  since  V(x)  >  0  for  all  points  in  R  ,  we  see  that  V(X(t))  must 

L/ 

converge  to  a  limit,  say  Cq,  as  t  goes  to  infinity. 

Let  L+  be  the  limiting  set  of  X(t) ,  that  is,  the  set  of  points 

p  in  E  for  which  there  exist  an  infinite  sequence  of  positive  times 
n 
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^l,^2’*’*,  such  that  ^(t^)  converges  to  p  as  k  goes  to  infinity. 
Note  that  the  limiting  set  must  be  an  invariant  set  of  (21)  .  For, 
suppose  L+  were  not  an  invariant  set.  Then  there  would  exist  a  solution 
X(t)  of  (21)  starting  for  t  =  0  at  a  point  in  L+  and  such  that 

X(t)  is  not  in  L+  for  t  >  0  and  t  small  enough.  Since  F(x)  is 
continuously  differentiable,  it  can  be  approximated  for  x  in  the  neigh¬ 
borhood  of  Pq  by  F(Pq)  +  K  •  (v  -  p  )  where  K  is  a  matrix  with 
constant  entries;  consequently,  for  x  near  enough  to  p^,  the  solutions 
of 

X  =  F(X) 

will  be  approximated  by  the  solutions  of 

X  =  F(Pq)  +  K  •  (X  -  pQ).  (23) 

However,  for  (23) >  it  is  easy  to  show  that  a  solution  such  as  X(t) 
cannot  exist  and  then,  by  continuity,  it  follows  that  X(t)  does  not 
exist  for.  (21).  Therefore,  L+  is  an  invariant  set. 

From  the  preceding  discussion  about  the  limit  of  V(X(t))  as  t 
goes  to  infinity,  it  follows  that  V(x)  =  for  all  x  in  L+; 

consequently,  V  =  0  on  L,.  Also,  the  set  L,  is  in  R_.  Therefore, 
if  M  is  the  largest  invariant  set  contained  in  the  set  of  points  at 
which  V  =  0,  then  M  contains  L  ,  the  limiting  set  of  the  solution 
X(t)  .  This  result  proves  the  theorem. 


It  is  useful  to  make  a  slight  generalization  of  Theorem  I  as  follows: 
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Theorem  II.  Suppose  for  all  x  in  E  that  V(x)  >  0  for 
x  f  0  and  that  V(x)  <  0;  suppose  further  that  V(x)  -*■  °°  as  x  -+  «>j 
then  every  solution  of  (21)  tends  to  the  set  M  of  Theorem  I. 

The  proof  follows  easily  from  Theorem  I.  Consider  any  solution 
X(t)  of  (21)  and  suppose  that  V(X(0))  =  C'  .  By  the  hypothesis  of  the 
theorem  the  set  of  points  in  such  that  V(x)  <  C  =  C'  +  £,  for 

arbitrary  e  >  0,  is  bounded.  Now,  we  can  use  Theorem  I  to  complete  the 
proof. 


Application  to  Networks 

From  equation  (20) ,  we  see  that  P(x)  will  be  a  Liapunov  function 
if  J  is  positive  definite  and  then  Theorem  II  could  be  applied  to  give  a 
stability  theorem.  In  general,  J  is  not  positive  definite  and  other 
arguments  must  be  used.  We  shall  conclude  with  a  statement  of  a  stability 
theorem  for  some  networks  [l]. 


Notice  that  if  is  a  matrix  and 

such  that  the  solutions  of  (19)  satisfy 


-  J  X  = 


9P1(x) 

9x 


P^(x)  a  scalar-valued  function 


(24) 


then,  because  of  linearity,  the  solutions  of  (19)  satisfy  all  equations 
of  the  form 


-  J  X  = 


_  9P(x) 
9x 


where 

J  =  aj  +  pj 
P  =  aP  +  pp 


(25) 
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with  a  and  fi  arbitrary  constants. 

A  possible  choice  for  J^,P^  ^4)  is  the  following: 

J1  =  PxxSJ’  P1  =  2  grad  P  ‘  S  grad  P>  C2&) 

where  S  is  a  symmetric  matrix  with  arbitrary  constants  for  entries. 

To  see  this,  start  from  (19)  ,  namely, 

-  J  X  =  grad  P  , 

and  multiply  in  front  by  the  matrix  P^S.  We  get 

-  J,  X  =  -  P  SJ  X  =  P  •  S  grad  P 

1  xx  xx  a 

but  notice  that 

|  g^Kgrad  P  •  S  grad  P)  =  |  P^  •  S  grad  P  +  ^  grad  P  •  SP^ 

=  Pxx  ’  S  grad  P ’ 

because  S  is  symmetric.  This  justifies  (26). 

Me  now  state  a  theorem  about  the  stability  of  networks.  For  the 
proof,  we  refer  to  [l,page  38]. 

Theorem  III.  Suppose 

P(x)  =  -  |(i,Ai)  +  B(v)  +  i  •  (yv  -  a), 

where  A  is  a  constant  symmetric  matrix,  y  is  an  arbitrary  constant 
matrix,  and  a  is  a  constant  vector.  Suppose  A  is  positive  definite. 
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B(  v)  +  |  YV  |  -*■  00  as  |v|  -*•  oo,  and  the  norm  of  the  matrix 

lJ(i)A-1-yC-ir(v) 

is  less  than  1-6,  with  6  >  0,  for  all  i,v;  then,  as  t  -*■  <»,  all 

the  solutions  of  (19)  tend  to  the  equilibrium  points,  that  is,  the  points 

at  which  grad  P  =0. 

The  idea  of  the  proof  is  to  find  a  suitable  matrix  S  in  ( 26) 

and  suitable  m  and  p  in  (25)  so  that  the  conditions  of  Theorem  II 

The  details  are  in  [l]. 


are  satisfied. 


AS 
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